Logarithmic boundary layers in highly turbulent Taylor-Couette flow 
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We provide direct measurements of the boundary layer properties in highly turbulent Taylor- 
Couette flow up to Ta = 6.2 x 10 12 using high-resolution particle image velocimetry (PIV). We 
find that the mean azimuthal velocity profile at the inner and outer cylinder can be fitted by the 
von Karman log law u + = — \ny + + B. The von Karman constant k is found to depend on the 
driving strength Ta and for large Ta asymptotically approaches n ~ 0.40. The variance profiles of 
the local azimuthal velocity have a universal peak around y + ~ 12 and collapse when rescaled with 
the driving velocity (and not with the friction velocity), displaying a log-dependence of y + as also 
found for channel and pipe flows Q], @| ■ 



Taylor-Couette (TC) flow is one of the paradigmatical 
flows in physics of fluids, next to Rayleigh-Benard con- 
vection, channel flow, and pipe flow. These flows have 
always been used in order to experimentally test new 
concepts in fluid dynamics. The TC system consists of 
two rotating coaxial cylinders shearing a fluid in between 
the cylinders, see fig. [TJ The Taylor number Ta quanti- 
fies the driving of this system and is directly related to 
the shear. For increasing Ta the system is first domi- 
nated by coherent structures Q whose length scale is of 
similar size as the gap width. For further increasing Ta 
turbulence develops in the bulk at length scales between 
the integral and the Kolmogorov scale while the bound- 
ary layers (BL) are still laminar. This regime, in which 
the flow has a turbulent bulk and the boundary layers 
are of Prandtl-Blasius type, is called the classical regime 
Q , and is analog to the situation in RB convection. By 
further increasing the Taylor number the system enters 
the so-called ultimate turbulent state in which also the 
boundary layers have turned turbulent These tur- 

bulent states are found to be remarkably similar for RB 
convection Q. Based on global transport measurements, 
the ultimate regime of turbulence sets in at Ra ~ 10 14 for 
RB H, but already at Ta - 10 s for TC [TM3|, due to 
the fact that mechanical driving in TC is much more ef- 
ficient in triggering BL instabilities than thermal driving 
is in RB. 

The coexistence of a laminar type boundary layer and 
turbulent bulk in classical turbulent RB convection has 
been well established from numerous experimental and 
numerical investigations and theoretical analysis 
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Very recently, a direct measurement [251 ] of 
the mean temperature profile close to the wall in the ulti- 
mate RB state revealed logarithmic behavior in the ulti- 
mate regime. Neither for TC nor for RB there had been 
any direct and systematical measurement of the velocity 
boundary layer in the highly turbulent ultimate state due 



to the experimental difficulties. However, the boundary 
layer properties are crucial to understand and define the 
general picture of the ultimate turbulent state. Assum- 
ing a logarithmic velocity profile in the boundary layers 
for highly turbulent TC flow, and matching the mean 
velocities at midgap, Lathrop et al. 

nam 

obtained a 

dependence of the global torque and the Reynolds num- 
ber, which agrees well with the torque measurements in 
ultimate turbulence regime 13, HI 23] ■ Only recently, di- 
rect measurements on BLs were conducted by van Hout 
& Katz [22| for Taylor number up to 2 x 10 9 , where they 
focused on the effect of counter-rotation and found that 
the von Karman constant depends on the angular veloc- 
ity ratio. 

In this letter we report the direct systematical exper- 
imental investigation of the boundary layer properties 
for very high Taylor numbers from Ta = 9.9 x 10 s to 
Ta = 6.2 x 10 12 using high resolution PTV and PIV 
2^-3^ with an unprecedented spatial resolution down 
to « 10 pm. We focus on the case of inner cylinder ro- 
tation, and examine the boundary layer properties as a 
function of Ta in the ultimate turbulent TC regime. 

In analogy to RB convection, in ref. [|| it was math- 
ematically found from the Navier-Stokes and continuity 
equations that the following quantity is strictly conserved 
in TC flow: 



Ju, = r 



(( u r^) z .e.t - vd r (^)z.e.t) ■ (!) 



Here {X) z g t t represents axial, azimuthal, and time aver- 
aging of X , u r is the radial velocity, oj the angular veloc- 
ity component uj = ug/r, and v is the kinematic viscosity. 
This flux is made dimensionlcss by dividing it by the flux 
for laminar flow: J^iam = 2vr 2 r 2 (uj i - u )/(r 2 - rf), 
giving an angular velocity Nusselt number: Nu w = 
Ju/ Ju,iam- This transport quantity is independent of r; 
any flux going through an imaginary cylinder with ra- 
dius r also goes through any other imaginary cylinder, or 
mathematically c^Nu^, = 0. This flux can be measured 
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FIG. 1: (color online) Sketch of the vertical cross section of the 
T 3 C [l]]]. The flow is illuminated from the side in the horizon- 
tal plane using a laser and the flow is imaged from the top us- 
ing a high-resolution camera. Top-right inset: Schematic top 
view of different regions inside the gap: IC (inner cylinder), 
OC (outer cylinder), and BL (boundary layer). Measurements 
are done at middle height. 



locally 27, 32| but also globally 1(J 33H35J by measur- 
ing the torque needed to sustain constant velocity of the 
cylinders. The torque T is related to the dimensionlcss 
torque G and to Nu w as follows: 



G 



T 



2ttI P v 2 



(2) 



where p is the density of the fluid, and £ the height of 
the cylinders. We can further relate these quantities to 
the wall shear stress t u , the friction velocity u T , and the 
viscous length scale 8 V at the inner cylinder wall: 



T 



(3) 

(4) 
(5) 



with an outer radius of i\ = 0.200 m, a transparent 
outer cylinder with inner radius r = 0.279 m, giving 
a radius ratio of 77 = 0.716. The cylinders have a 
height of I = 0.927 m, resulting in an aspect ratio of 
T = £/(r — Ti) = 11.7. More details can be found in 
ref. [3l|. In this letter the focus is on inner cylinder ro- 
tation only, the outer cylinder remains fixed during all 
the experiments. For several rotation rates we perform 
PIV and PTV (particle tracking velocimetry) measure- 
ments in our setup. The working fluid (water) is seeded 
with fluorescent polymer particles [36| with a diameter of 
1-20 |J.m. Using a laser [37[ we create a horizontal laser 
sheet of roughly 500 \mi thickness that illuminates the 
particles inside the fluid. The flow is then imaged from 
the top, see also figure [TJ using a high- resolution camera 
[38| with large dynamical range. In front of the camera a 
50 mm objective lens was mounted, with the working dis- 
tance of approximately 600 mm, giving a scaling factor 
of ?s 54 |xm/px. 
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FIG. 2: (color online) Azimuthal velocity profiles for varying 
Ta across the gap of the Taylor-Couette apparatus. The leg- 
end indicates the Taylor numbers of the experiments. Insets a 
and b show a zoom of the data of the inner and outer bound- 
ary layer, respectively. Individual data-points are plotted in 
all 3 figures, showing the high resolution of the measurements. 
The gray solid line is the exact laminar solution of the Navier- 
Stokes equations for infinitely long cylinders. 



Note that as Nu^ is conserved radially, it is the same 
at both cylinders, and using eq. ([2]), we see that thus 
also the torque T at both cylinders should be the same. 
Consequently, t u , u Tl and <5„ are different at the inner 
and outer cylinder, and the following relations hold: 



tw/t^o = 1/rf- 
u T ,i/u T ^ = l/l] 
Su,i/Sv,o = V 



(6) 

(8) 



The apparatus used for the experiments, the Twente 
turbulent Taylor-Couette (T 3 C), has an inner cylinder 



For each rotational velocity 10 image pairs were ac- 
quired at a recording frequency of 10 Hz. The mean 
velocity distribution was computed using single-pixel 
ensemble-correlation, leading to a final resolution of 
w 150 |_im, resulting in over 500 independent datapoints 



in the 80 mm gap [28|, |39( . The standard deviation was 
directly computed from the velocity probability density 
function, which was extracted from the shape of the cor- 
relation function, as discussed in ref. [3(|. This proce- 
dure ensures that all turbulent scales are included in the 
standard deviation. In contrast to standard PIV analysis, 
where only spatially low-pass filtered results arc achieved, 
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FIG. 3: (color online)(a) Azimuthal velocity profile near the inner cylinder (r' G [0, 1/2]) for varying Ta. y + is the distance 
from the inner cylinder in units of the viscous length scale 5 V . u + is defined as (u(r'j) — u(r)) /u T , where w(r) is the azimuthal 
component of the velocity and u{ri) is the azimuthal velocity of the inner cylinder, (b) Azimuthal velocity profile near the 
outer cylinder r' £ [1/2, 1]. y + = (r — r )/S v is the distance from the outer cylinder in wall units. For the case of the outer 
boundary the velocity is scaled as u + = u(r)/u T . Both figures include the logarithmic law of the wall u + = - In y + + B by von 
Karman, with the typical values of k = 0.40 and B = 5.2, the viscous boundary layer u + = y + , and in gray the fitting domain 
y + £ [50, 600] . Each data-set is accompanied by an arrow with caption indicating the Taylor number. 



here also the contribution of the small scale fluctuations 
are properly considered. In order to resolve the near wall- 
region at the inner cylinder, a microscope (40| with a 
focal length of 300 mm was mounted in front of the cam- 
era. With this setup, a scaling factor of « 10 um/px was 
achieved. We evaluate the near-wall region with PTV 
methods, which is well suited for this purpose (29[. 

Figure [2] shows all the measured profiles; 5 covering 
the full gap, and 2 covering just the region near the in- 
ner cylinder. We see that the profiles do not conform 
to the laminar solution (for T = oo) and that the bulk 
has a much shallower slope — as seen in RB convection 
[i~5| — indicating turbulent mixing. For the laminar solu- 
tion the convective part of is zero; any convection in 
the bulk of the system therefore decreases the d r oj term, 
resulting in a shallower angular velocity gradient. From 
the insets we can clearly see that for increasing Taylor 
number the boundary layers become steeper and steeper, 
and indeed the angular velocity profile has then to be- 
come steeper as u r of eq. ((J) is zero at the wall, and we 
are only left with the term d r u>. We now split our data 
in two parts: r' G [0,1/2] (inner BL) and r' € [1/2,1] 
(outer BL) and normalize velocities with the appropriate 
u T (eqs. ((4|) and ([7])) and distances with the respective 
b v (eqs. ([5]) and (Qp), see figure [3] We used global torque 
measurements [111 . |34| ] to find u T and 8 V . In the viscous 
sublayer (y + < 5) the velocity profiles are u + = y + which 
we can nicely resolve. It should be noted that global 
torque measurements provide an average torque r, while 
the local torque depends on height (following the large- 
scale Taylor vortex structure) . So while eq. (j3J) holds for 



the average r w .i, it might have an axial dependence. To 
normalize the velocity profiles we used the globally mea- 
sured torque and therefore we have the average of t w .;, 
this causes the imperfect matching of u + = y + . The mea- 
surements were performed at middle height as sketched 
in fig. □ 

For y + > 50 the effects of viscosity diminish — indeed 
the contribution of vd r ui to J w is very small in the bulk of 
the flow. Furthermore, as suggested by Prandtl and von 
Karman 4l|, in this limit the velocity profile converges 
to: 



,+ - 



1 



lny+ + B, 



(9) 



with k the von Karman constant and B the logarithmic 
intercept. 

Indeed, outside the viscous wall region y + > 50 (inside 
the outer layer) figure [3] conclusively shows the existence 
of a log-layer in the ultimate TC regime (Ta > 10 9 ), 
which is in sharp contrast to the laminar boundary lay- 
ers found in the DNS simulations 23j (indicated by thin 
lines) in the classical turbulent regime. Our profiles are 
fitted to eq. ([9]) over the range y + € [50, 600] (see the 
shaded area in figurc[3|) and the values for k are extracted, 
sec figured] They are found to be different for the inner 
and outer BL and depend on Ta. For increasing Taylor 
number we see that the values for the inner BL approach 
a value of k « 0.40, close to the known classical value 
of k = 0.40 @, S3], recently systematically examined by 
Marusic et al. [l| . We included in figure |4] a line indi- 
cating the influence of the curvature in the log-regime 
(600(5„/r). It can be observed that for our lower Taylor 
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FIG. 4: (color online) The parameter re obtained by fitting 
u + = - In y + + B for y + £ [50, 600] near the inner and outer 
cylinder. The color scheme is identical to figures [2] and [3] the 
black symbols are DNS results [2^1 . The solid colored line 
shows the ratio of the extent of the log-layer and the radius 
of the cylinder, indicating the influence of the curvature. 



numbers the importance of the curvature cannot be ne- 
glected, and that only for the very high Taylor numbers 
the curvature becomes negligible (green part of the solid 
line). For Ta > 10 13 we expect the effect of curvature 
to become negligible and that the value for k approaches 
the classical value [l|. We also would like to note that 
due to possible height dependence of t u ^, and therefore 
the imperfect matching of u + = y + , the values of k could 
slightly vary. 

In addition to PIV, also high resolution PTV measure- 
ments and analysis have been performed. For these mea- 
surement we zoomed into the area near the inner cylinder 
using a long-distance microscope to obtain a scaling fac- 
tor of w 10 pm/px. The spatial resolution of PTV only 
depends on the number of images, and can thus be better 
than the pixel grid spacing projected in to physical space 
[28j| . We extract the variance a 2 (ug) from the shape of 
the probability density function of the correlation func- 
tion [3(|, see figure El We normalize a(ue) 2 with the 
driving velocity u(r.;) (see fig. [5^,) and with the friction 
velocity u T (see fig. [SJd). For both curves the maximum 
of a(ug) 2 is around y + = 12, which is remarkably similar 
to the values obtained in pipe and channel flows (see e.g. 
[42L I43I ] ) . These peaks universally collapse inside the vis- 
cous wall region (y + < 50), where curvature effects do not 
play a role like in the log-regime. In addition it can be ob- 
served that indeed the collapse of the data is better when 
we normalize a(ug) with the driving velocity rather than 
the shear velocity u T . As shown in fig.OJi, we fit the data 
for y+ > 100 with (a(u g )/u T ) 2 = B x - A 1 \ny+ —the 
log law for the velocity variance 0. The corresponding 
fitting parameters are indicated in fig. [SJi. Remarkably, 
the slope A\ (varying from 1.22 to 1.53) is comparable 
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FIG. 5: (color online) (a) The variance of the local azimuthal 
velocity is presented as a function of radial distance from the 
inner cylinder. The velocity is made dimensionless using the 
friction velocity u T , and the distance is in wall-units y + = 
(r — n)/6u. (b) Same as in figure a but the velocity has 
been rescaled using the driving velocity u(ri). Corresponding 
colors in figures a and b correspond to the same data-set, 
and are consistent with previous figures. The Taylor numbers 
are indicated with arrows in figure a. For y + > 100 the 4 
highest Ta cases are fitted with (a(ug)/u T ) 2 = Bi — Ai lnj/ + . 
The fitting parameters are listed inside figure a and colored 
accordingly. 



with values found (^4i~ 1-25) in high Reynolds number 
boundary layer flows [2|. 

In conclusion, wc performed direct measurements of 
the velocity boundary layer profiles in highly turbulent 
Taylor Couette flow up to Ta = 6.2 x 10 12 . In contrast 
to the laminar boundary layers in the classical turbulent 
TC regime, the present data in the ultimate regime pro- 
vides direct experimental evidence of the emergence of a 
log-layer, as theoretically proposed in refs. 0,0- The 
fitted von Karman constant k is found to approach the 
classical value of 0.40 for large enough Ta. Furthermore, 
we find that the peak in the standard deviation of the 
azimuthal velocity universally collapses around y + = 12, 
and that the height of the peak is found to collapse bet- 
ter when scaled with the driving velocity as compared to 
the friction velocity. Lastly, the variance profiles depict 
a log dependence for y + > 100. 
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